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Numerical simulation results are presented which suggest that a class of non-adiabatic 
rapid passage sweeps first realized experimentally in 1991 should be capable of imple- 
menting a universal set of quantum gates Q u that operate with high fidelity. The gates 
constituting Q u are the Hadamard and NOT gates, together with variants of the phase, 
7r/8, and controlled-phase gates. The universality of Q u is established by showing that it 
can construct the universal set consisting of Hadamard, phase, 7r/8, and controlled-NOT 
gates. Sweep parameter values are provided which simulations indicate will produce the 
different gates in Q u , and for which the gate error probability P e satisfies: (i) P e < 10~ 4 
for the one-qubit gates; and (ii) P e < 1.27 X 10 — 3 for the modified controlled-phase gate. 
The sweeps in this class are non-composite and generate controllable quantum interfer- 
ence effects that allow the gates in Q u to operate non-adiabatically while maintaining 
high fidelity. These interference effects have been observed using NMR, and it has pre- 
viously been shown how these rapid passage sweeps can be applied to atomic systems 
using electric fields. Here we show how these sweeps can be applied to both supercon- 
ducting charge and flux qubit systems. The simulations suggest that the universal set 
of gates Q u produced by these rapid passage sweeps shows promise as possible elements 
of a fault-tolerant scheme for quantum computing. 

Keywords: fault-tolerant quantum computation, accuracy threshold, quantum interfer- 
ence, resonance, non-adiabatic dynamics, superconducting qubits 

1 Introduction 

Slightly more than a decade has passed since the accuracy threshold theorem was first proved, 
establishing the unexpected result that, under appropriate conditions, a quantum computation 
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of arbitrary duration could be carried out with arbitrarily small error probability in the 
presence of noise, and using imperfect quantum gates pQ-[E|- The required conditions are 
that: (1) the computational data is protected by a sufficiently layered concatenated quantum 
error correcting code; (2) fault-tolerant protocols for quantum computation, error correction, 
and measurement are used; and (3) all quantum gates used in the computation have error 
probabilities (per operation) P e that fall below a value known as the accuracy threshold 
P a . This threshold has been calculated for a number of simple error models with results 
falling in the range 1CP 6 < P a < 1CP 3 . For many, P a ~ 10~ 4 has become a rough-and- 
ready working estimate for the threshold so that gates are anticipated to be approaching 
the accuracies needed for fault-tolerant quantum computing when P e < 10~ 4 . A number 
of universal sets of quantum gates have been found [5]-[]3]. With such a set of gates, any 
./V-qubit unitary operation can be applied using a quantum circuit made up solely of gates 
belonging to the universal set. Thus the problem of producing sufficiently accurate quantum 
gates shifts to producing a sufficiently accurate universal set of such gates. One well-known 
universal set consists of the one-qubit Hadamard, phase, and tt/8 gates, together with the 
two-qubit controlled-NOT gate [14]. 

In this paper numerical simulation results are presented which suggest that a class of 
non-adiabatic rapid passage sweeps known as twisted rapid passage (TRP), first realized 
experimentally in 1991 [15], should be capable of implementing a universal set of quantum 
gates Q u that operate non-adiabatically and with high fidelity. Q u consists of the one-qubit 
Hadamard and NOT gates, together with variants of the: (i) one-qubit phase and 7r/8 gates; 
and (ii) two-qubit controlled-phase gate. The universality of this set is established by showing 
that it can construct the universal set of Ref. |14j . Sweep parameter values are provided 
which simulations indicate will produce the different gates in Q u , and for which the gate error 
probability P e satisfies: (i) P e < 1CT 4 for the one-qubit gates; and (ii) P e < 1.27 x 10~ 3 for 
the modified controlled-phase gate. This level of accuracy is a consequence of controllable 
quantum interference effects that are generated by this class of rapid passage sweeps [Ml [T7l 
118] . The sweep parameter values provided for each gate were found using an optimization 
procedure that searches for minima of P e . TRP implementation of the one-qubit gates in Q u 
was shown in Ref. [18]. Here we complete the universal set Q u by showing that TRP can 
also implement a high fidelity, non-adiabatic modified controlled-phase gate. TRP sweeps 
have already been experimentally realized in NMR systems [H] [T7] . Ref. [18] showed how 
TRP sweeps can be applied to atomic systems using electric fields. Here we extend the list of 
physical systems to which TRP can be applied by showing how these sweeps can be applied 
to both superconducting charge and flux qubit systems. 

The outline of this paper is as follows. Section [2] begins with a summary of the essential 
properties of TRP. It then details our simulation protocol, and the procedure used to optimize 
the TRP sweep parameters. Section [3] presents our simulation results for the universal set 
of gates Gu- For each gate in Q u , we present: the optimized sweep parameter values that 
produce it; an upper bound on the gate error probability P e ; and the gate fidelity T. For the 
two-qubit modified controlled-phase gate, we also include the unitary operation produced. 
We have not included the corresponding unitary operations for each of the one-qubit gates in 
Q u as they appear in Ref. [TS] . Section 2] then explains how TRP sweeps can be applied to 
both superconducting charge and flux qubits. Finally, Section [5] makes closing remarks. 
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2 Preliminaries 

This section discusses a number of important background topics. Section [2.11 provides a brief 
review of TRP and shows how these sweeps generate quantum interference effects that are 
controllable through variation of the sweep parameters. Section 12.21 describes the protocol 
used to simulate the one- and two-qubit gates that are of interest in this paper. Finally, to 
produce high performance gates, the simulations are combined with an optimization procedure 
described in Section 12.31 that searches for sweep parameter values that minimize an upper 
bound on the gate error probability P e . 

2.1 Twisted Rapid Passage 

To introduce TRP we consider a single qubit interacting with an external control field F(£) 
via the Zeeman interaction 

H z (t) = -*-F(t), (1) 

where er, are the Pauli matrices. The sweeps we will be interested in are a generalization of 
those used in adiabatic rapid passage (ARP) [19] . In ARP the field F(t) in the detector frame 
[2"0] is inverted over a time To such that F(i) = bx + atz. In the NMR realization of ARP, as 
seen in the lab frame, the detector frame rotates about the static magnetic field Bq z. In the 
detector frame, z is chosen parallel to the rotation axis. In the rotating wave approximation 
the rf-magnetic field B r j in the lab frame lies in the x-y plane and rotates about the static 
magnetic field. The detector frame is chosen to rotate with B r y so that, in this frame, the rf- 
field is static and its direction defines x: B^f = bit. The inversion time To is large compared 
to the inverse Larmor frequency uJq 1 (viz. adiabatic), though small compared to the thermal 
relaxation time Ttn (viz. rapid). ARP sweeps provide a highly precise method for inverting 
the Bloch vector Sj = (<Tj), although the price paid for this precision is an adiabatic inversion 
rate. We are interested in a type of rapid passage in which the control field F(t), as seen in 
the detector frame, is allowed to twist around in the x-y plane with time-varying azimuthal 
angle (f>(t), while simultaneously undergoing inversion along the z-axis: 

F(i) = b cos(j)(t)x + b sin0(<)y + atz . (2) 

Here —To/2 < t < Tq/2 and y = z x x. Note that any pair of orthogonal unit vectors in 
the x-y plane can be used for x and y. Different choices simply alter the value of <fi(t = 
0). As will be seen shortly, interesting physical effects arise when the twist profile cj)(t) 
is chosen appropriately. This class of rapid passage sweeps is referred to as twisted rapid 
passage (TRP). The first experimental realization of TRP in 1991 by Zwanziger et al. [15] 
carried out the inversion adiabatically with <p(t) = Bt 2 . Subsequently, non- adiabatic TRP 
was studied with polynomial twist profile 4>(t) = (2/n)Bt n [IB], and controllable quantum 
interference effects were found to arise for n > 3. Zwanziger et al. [17j implemented non- 
adiabatic polynomial TRP with n — 3,4 and observed the predicted interference effects. In the 
Zwanziger experiments [1 51 117] , a TRP sweep is produced by sweeping the detector frequency 
4>det{t) linearly through resonance at the Larmor frequency loq: 4>det(t) = ujq + (2at)/h. The 
frequency of the rf- field 4> r f(t) is also swept through resonance in such a way that 4> r f{t) = 
4>det{t) — 4>(t), where <j>(i) = (2/n)Bt n is the TRP polynomial twist profile. Substituting the 
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expression for 4>det(t) into that for 4> r f(t) gives 

^r/W = wo + y- #) ■ (3) 
At resonance 4>rf(t) = i^o- Inserting this condition into eq. ((3]), it follows that at resonance 

at-~0(t)=O . (4) 

For polynomial twist <fi(t) = (2/n)Bt n , eq. J3|) has n — 1 roots, though only real-valued roots 
correspond to resonance. Ref. 16J showed that for n > 3, multiple passes through resonance 
occur during a single TRP sweep: (i) for all n when B > 0; and (ii) for n odd when B < 0. 
We restrict ourselves to B > in the remainder of this paper. In this case, the qubit passes 
through resonance at the times: 

f 0, (a/UB)^ (n odd) , . 

\ 0,±(a/hB)^ (n even) 

We see that the time separating the qubit resonances can be altered by variation of the sweep 
parameters B and a. Ref. [16] showed that these multiple resonances have a strong influence 
on the qubit transition probability. It was shown that qubit transitions could be significantly 
enhanced or suppressed by small variation of the sweep parameters, and hence of the time 
separating the resonances. Plots of the transition probability versus time suggested that the 
multiple resonances were producing quantum interference effects that could be controlled by 
variation of the TRP sweep parameters. In Ref. [21j the qubit transition amplitude was 
calculated to all orders in the non-adiabatic coupling. The result found there can be re- 
expressed as the following diagrammatic series: 

Lower (upper) lines correspond to propagation in the negative (positive) energy level and the 
vertical lines correspond to transitions between the two energy levels. The calculation sums 
the probability amplitudes for all interfering alternatives |22j that allow the qubit to end up in 
the positive energy level at time t given that it was initially in the negative energy level. As we 
have seen, varying the TRP sweep parameters varies the time separating the resonances. This 
in turn changes the value of each diagram in eq. ([6]), and thus alters the interference between 
alternatives in the quantum superposition. Similar diagrammatic series can be worked out for 
the remaining 3 combinations of final and initial states. It is the sensitivity of the individual 
alternatives/diagrams to the time separation of the resonances that allow TRP to manipulate 
this quantum interference. Zwanziger et al. [17j observed these interference effects in the 
transition probability using liquid state NMR and found quantitative agreement between 
theory and experiment. It is the link between the TRP sweep parameters and this quantum 
interference that we believe makes it possible for TRP to drive highly accurate non-adiabatic 
one- and two-qubit gates. The results presented in Section [3] for the different gates in the 
universal set Q u are found by numerical simulation of the Schrodinger equation. We next 
describe how the simulations are done. 
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2.2 Simulation Protocol 

As is well-known, an iV-qubit quantum gate applies a fixed unitary transformation U to 
iV-qubit states: 

\^out) = U\ip m ). (7) 

We will be interested in the unitary transformations applied by the: (i) one-qubit Hadamard 
(U H ), phase (Up), tt/8 (U„ /6 ), and NOT (U N OT ) gates 
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and (ii) the two-qubit controlled-NOT (U< 
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-) and controlled-phase (U CP ) gates 
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The one- and two-qubit matrices appearing in eqs. ([5|)- (|10p are in the representation spanned 
by the one- and two-qubit computational basis states |i) and which are, respectively, the 
eigenstates of a z and a\ ® <j\ : 



= (-i) 



1 °i 



\ij) (i, j = 0,l). 



(11) 



The dynamical impact of TRP is determined by numerical simulation of the Schrodinger 
equation. Let H(t) denote the Hamiltonian for an iV-qubit system, and \Ek(t)) the instan- 
taneous energy eigenstates which satisfy H(t) \Ef.(t)) = Ek(t) \Ef.(t)) with k = 1, . . . , 2^. It 
is found that the numerical stability of the simulations is enhanced if we expand \ip(t)) in 
the instantaneous energy eigenstates \Ek(t)). Because of the direct connection between these 
states and H[t), they carry substantial dynamical information, and a substantial portion of 
the dynamics due to H(t) can be accounted for by choosing this basis. This makes the task 
of determining the remaining dynamics using the Schrodinger equation much simpler and the 
simulations more stable. We thus write 



\iP(t)) = a k(t) \E k (t)) exp 



k=l 



T /2 



dT{E k {T)-hj k (T)} 



(12) 



where 7fc(£) is the adiabatic geometric phase [23] associated with the energy level E k (t), and 



i{E k (t)\-\E k (t)). 



(13) 



Substituting eq. (|12[) into the Schrodinger equation, and using the orthonormality of the 
instantaneous energy eigenstates, one arrives at the equations of motion for the expansion 
coefficients a k (t): 



-|r = ~Y1 a K*) r w(*) ex P 

Z^fc 



To/2 



dr 5i k (r) 



(14) 
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where k = 1, . . . , 2 , and 



r H (t) = (B fc (t)|-|£,(t)) 



(15) 



M*) = ^I^ -CtK*)-^)). (16) 



The simulations set the initial AT-qubit state to be one of the computational basis states 
\ix ■ ■ ■ i 2 n). The simulation outcome is the \i± ■ ■ ■ i 2 w)-column of the unitary transformation U 
produced by the TRP sweep. By simulating all 2 N computational basis states, U is determined 
column- by-column. 

It proves useful to recast eqs. (| 14[) in dimensionless form. For the one-qubit case we 
introduce the dimensionless time r = (a/b)t, the dimensionless inversion rate A = ha/b 2 , and 
the dimensionless twist strength rj n = (fiB /a)(b/a) n ~ 2 . For the remainder of this paper we 
restrict ourselves to quartic TRP which has twist profile <fi n (T) with n = 4: 



Mr) = \B 

2 V a 



(17) 



The sweep parameters a, 6, and B are chosen to be positive so that A and r; 4 are also positive. 
The Hamiltonian for the one-qubit gates is the Zeeman Hamiltonian Hz(t) given in eqs. ([T|) 
and ([!]). The dimensionless version of Hz(t) is found by multiplying it by (b/fia). Denoting 
the result by Hi(t), one finds 



r 1 

-Hi(t) = "T^ - J [cos^cr^ 



sin <p4 a y 



(18) 



where ^(r) is given by eq. (|17[) . H\(t) is the Hamiltonian that drives the one-qubit simula- 
tions via eqs. (| 14[) (p~6j) in dimensionless form. 

The derivation of the dimensionless Hamiltonian H2(t) that drives the two-qubit sim- 
ulations is more complicated and is presented in Appendix A. As explained there, H^ir) 
includes an interaction that Zeeman-couples each qubit to a TRP control field, as well as 
an Ising interaction between the two qubits. The Ising interaction was chosen because of 
its simplicity, and because of its occurrence in many physical systems. Modification of the 
simulations to include alternative qubit-qubit interactions is straightforward. The resulting 
two-qubit Hamiltonian yields energy level spacings that give rise to a degeneracy in the reso- 
nance energy for the energy level pairs (E\ <-> E2) and (E3 <-> E4). To remove this degeneracy 
the term AH = C4 \E4,(t)) (E4,(t)\ is added to the Hamiltonian, where C4 is a constant. The 
end result is the Hamiltonian H%{t) that drives the two-qubit simulations: 



H 2 (t) 



(di+d 2 ) t 
2 A 

(h r 
2 A 



a\ z - y [ cos 4>i a lx + sin 4 a ly ' 



&2z 



A 



COS 04 2x + SUl 04 CJ 2 y 



-d 4 u\ z a 2z + AH. 



(19) 



The constants di (i — 1, . . . , 4) are defined in Appendix A where a description of their physical 
significance is also given. 
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Having described how the simulations are done, and how we determine the actual unitary 
transformation U a produced by a specific assignment of the TRP sweep parameters, we go 
on in Section [2. 31 to explain how the sweep parameters are iteratively modified so as to make 
U a approach a target gate U t as closely as possible. 

2.3 Sweep Parameter Optimization 

Let 7ijv be the Hilbert space for an 7V-qubit system. As in Section 12. 21 let U a denote the 
actual unitary operation produced by a given set of TRP sweep parameters, and Ut a target 
unitary operation we would like TRP to approximate as closely as possible. Introducing the 
operators D = U a — Ut and P = D^D, and the normalized state \tp), we define \ijj a ) — U a \tp), 
\if)t) = U t \4>), and \ip±) = {I — |V't)('0*| HV'a)- Ref- [IB] showed that the error probability 
P e (V>) for U a acting on the state \tp) is Pe{ip) — ±_\^p ±)\ 2 . The gate error probability P e for 
the TRP gate U a was then defined to be the worst-case error probability: 

P e = max P e {i>). (20) 

W 

Because of the search over evaluation of P e from eq. (f2T))) is not practical. Ref. [H] showed 
that 

P e <TrP : (21) 

where P = D^D was defined above. Once U a has been determined by the numerical sim- 
ulation, Tr P is easily calculated, and so it makes a convenient proxy for P e in the sweep 
parameter optimization procedure that we are now ready to describe. 

To find TRP sweep parameter values that yield highly accurate non-adiabatic quantum 
gates it proved necessary to combine our simulations with function minimization algorithms 
[24] that search for sweep parameters that minimize the upper bound Tr P for the gate error 
probability P e . The multi-dimensional downhill simplex method was used for the one-qubit 
gates in the universal set Q Ul while simulated annealing was used for the two-qubit modified 
controllcd-phase gate. For quartic TRP, the sweep parameters are (A, 774) which can be 
thought of as specifying a point in a two-dimensional parameter space. The downhill simplex 
method takes as input three sets of sweep parameters which specify the vertices of a simplex 
in the two-dimensional parameter space. The dynamical effect of the TRP sweep associated 
with each vertex is found by numerically integrating the Schrodinger equation as described in 
Section [2. 21 The output of the integration is the unitary operation U a that a particular sweep 
applies. Having U a we determine P — (U a — Uty (U a — Ut) and evaluate Tr P. The downhill 
simplex method iteratively alters the simplex (i.e. one or more of its vertices) until sweep 
parameter values are found that yield a local minimum of Tr P. Because the minimum found 
is not global, some starting simplexes will give deeper minima than others. Though there was 
no guarantee, it was hoped that a starting simplex could be found that yielded Tr P < 10 -4 . 
Some trial and error in specifying the starting simplex was thus required, though for the one- 
qubit gates in Q u , the trial and error procedure eventually proved successful. For the two-qubit 
modified controlled-phase gate, trial and error proved inadequate, and it was necessary to use 
simulated annealing to find sweep parameter values that yielded Tr P ~ 10 -3 . We present 
our results for all gates in Q u in the following section. 
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3 Universal Set of Quantum Gates 

To determine fidelities for the TRP gates we use the fidelity 



F n = -Be[Tr(UlU t )], 



(22) 



where n denotes the number of qubits acted on by the gate. T n is an extension to n qubits 
of the fidelity used in Ref. [25 for one-qubit gates (n = 1). It is possible to relate our Tr P 
upper bound for P e to T n . Recalling that P — (U a — U t y (U a ~ U t ), it follows that 



Tr P = Tr 



- (ulut + u}u a ) 

2 n+i _ 2 Re [Tr (U%U t )] 
2 n+1 (1 - T n ) , 



and so 



Tn = 1 - 



1 



2 n+l 



Tr P. 



Tr P thus yields the fidelity T n and an upper bound on the gate error probability P e . 
Next, notice that the gates Up and J7 T /g (see eqs. ([5]) and ©) can be re-written as 



where 



Up 

Utt/8 



v P 



e"/ 4 U NOT V P 
e n/8 U NOT V„ /8 , 





i7r/4 



e 




iir/A 



i7r/8 







(23) 
(24) 



(25) 
(26) 



(27) 
(28) 



and U NOT is given in eq. ©. Note also that U CNOT (eq. (jTUJ)) can be re-written as 

U CNOT = (I 1 ®U%)[( a\ ® I 2 ) V CP ] (I 1 ® U 2 H ) . (29) 

Here the superscript on a one-qubit gate labels the qubit on which the gate acts, and we have 
introduced the modified controlled-phase gate 



V a 



(1 

1 



\ 



\ 



-1 

1 / 



(30) 



Finally, notice that the Pauli matrix <j l z can be implemented using the phase gate U l P : o\ = 
(U P ) 2 . It follows from this and eqs. (US]), P5|l. and that the set of gates {U H , U P , U v/8 , 
U C not } can be constructed using the set Q u = {{/#, U NOT , Vp, V^/g, V CP }. Since the first set 
of gates is universal [TJ], so is the set Q u . As will be seen below, TRP can be used to produce 
all gates in Q u . For each gate we present our best-case result and show how gate performance 
is altered by small variation of the parameters. 
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3.1 One-Qubit Gates 

A study of the TRP-implementation of the one-qubit gates in Q u was first reported in Ref . [18] . 
The essential results are included here for the reader's convenience, though for the sake of 
brevity, we have not reproduced the unitary operation U a generated by TRP for each gate. 
The interested reader can find them displayed in Ref. [Tg] . 

As noted in Section 12. 2\ all results presented in this paper are for quartic TRP which has 
twist profile 

Here r, A, and r\ t are dimensionless versions of the time t, inversion rate a, and twist strength 
B. For the reader's convenience, their definitions (Section l2.ip are restated here: 



(- 
Kb 



ha ( hb 2 " , , 

7^ 5 '/. (—)/'• IV) 



Throughout this paper we assume a, 6, and B are positive. The parameter b was introduced 
in eq. ^ and is proportional to the rf-field amplitude in an NMR realization of TRP [TBI HZ] ■ 
All simulations were done with A > 1 corresponding to non-adiabatic inversion [TH] [T7j , and 
the dimensionless inversion time tq = aT^/b was fixed at 80.000 for the one-qubit simulations, 
and at 120.00 for the two-qubit simulations. 

The translation key connecting the simulation parameters to the experimental sweep pa- 
rameters used in the (one-qubit) Zwanziger experiments [151 117j is given in the Appendix 
of Ref. [16]. We re- write the formulas for quartic twist here for convenience. Note that 
Zwanziger's symbol B is here replaced by B to avoid confusion with our use of B to denote 
the twist strength. First we give the formulas connecting our parameters (a, b, B, To) to the 
Zwanziger parameters (lji, A,B,Tq): 

ui = T (33) 
n 

A = (34) 

B = (35) 

where the inversion time Tq is common to both parameter sets. The formulas linking the 
dimensionless sweep parameters (A, r\A to the Zwanziger parameters (lui,A, B,T ) are: 



AA 
Buj 2 

2A 3 Tn 



A = -277T (36) 



(37) 



In the experiments of Ref. [17]: = 393 Hz; T = 41.00 ms; A = 50 000 Hz; and B was 
calculated from eq. ([57]) with rj 4 varying over the range [4.50,4.70] x 10~ 4 . 

Hadamard Gate 

Here the target gate Ut is the Hadamard gate Uh (see eq. ©). The TRP sweep parameters 
A = 5.8511 and tj 4 = 2.9280 x 10~ 4 produce a unitary gate U a [TS] for which Tr P = 8.82 xl0~ 6 . 
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This yields a gate fidelity T H = 0.9999 98 and the bound P e < 8.82 x 10" 6 on the gate error 
probability. Table [1] shows how gate performance varies when the sweep parameters are 
altered slightly. Of the two sweep parameters, t] 4 is seen to have the largest impact on gate 



Table 1 . Variation of Tr P for the Hadamard gate when the TRP sweep parameters are altered 
slightly from their best performance values. The columns to the left of center have 774 = 2.9280 X 
10 — 4 and those to the right have A = 5.8511. 





A 


TrP 


A 






TrP 




2.9280 x 10~ 4 


5.8510 


7.22 x 10- 5 


5.8511 


2.9279 x 10" 


-4 


7.03 x 10 


-4 




5.8511 


8.82 x 10~ 6 




2.9280 x 10" 


-4 


8.82 x 10 


-6 




5.8512 


1.84 x 10" 5 




2.9281 x 10" 


-4 


6.14 x 10 


-4 



performance. This will turn out to be true for the other one-qubit gates as well. Although 
TRP can produce a non-adiabatic Hadamard gate whose error probability falls below the 
accuracy threshold P a ~ 10~ 4 , it is clear from Table Q] that the sweep parameters must be 
controlled to 5 significant figures to achieve this level of performance. See Section[5]for further 
discussion. 

V P Gate 

The modified phase gate Vp (see eq. (f27| ) is the target gate here. The sweep parameters A = 
5.9750 and r? 4 = 3.8060 x 10~ 4 produce a unitary gate U a [TSj which has TrP = 8.20 x 10~ 5 . 
This gives a fidelity T Vp = 0.9999 80 and the bound P e < 8.20 x 10~ 5 . Table EH shows how 
TrP varies with small changes in A and r\ A , Again, gate performance is most sensitive to 



Table 2. Variation of TrP for the modified phase gate Vp when the TRP sweep parameters 
are altered slightly from their best performance values. The columns to the left of center have 
774 = 3.8060 X 10~ 4 and those to the right have A = 5.9750. 



V* 


A 


TrP 




A 






TrP 




3.8060 x 10" 4 


5.9749 


1.56 x 10" 


4 


5.9750 


3.8059 x 10" 


-4 


2.29 x 10 


-3 




5.9750 


8.20 x 10" 


5 




3.8060 x 10- 


-4 


8.20 x 10 


-5 




5.9751 


1.43 x 10" 


4 




3.8061 x 10" 


-4 


1.88 x 10' 


-3 



variation of 774, and the sweep parameters must be controlled to high precision to surpass the 
accuracy threshold (see Section [SJ . 

Vx/8 Gate 

The target gate this time is the modified ir/8 gate Vwg (see eq. (1281) ). For A = 6.0150 and 
i] 4 = 8.1464 x 10" 4 , TRP produces a unitary gate U a [18] which has TrP = 3.03 x 10~ 5 , 
fidelity Tv^ /e , = 0.9999 92, and P e < 3.03 x 10 -5 . Table [3] shows how gate performance varies 
when the sweep parameters are altered slightly. As with the two previous gates, performance 
is most sensitive to variation of r) 4 , and the sweep parameters must be controllable to high 
precision (see Section [5]). 

NOT Gate 

Here the target gate is the NOT gate (see eq.©). For A = 7.3205 and 7/4 = 2.9277x 10~ 4 , TRP 
produces the unitary gate U a [18j for which TrP = 1.10 x 10~ 5 , fidelity T NO t = 0.9999 97, 
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Table 3. Variation of Tr P for the modified it/8 gate when the TRP sweep parameters are altered 
slightly from their best performance values. The columns to the left of center have r\& = 8.1464 X 
10~ 4 and those to the right have A = 6.0150. 



V* 


A 


TrP 




A 






TrP 




8.1464 x 1CT 4 


6.0149 


1.30 x 10" 


:i 


6.0150 


8.1463 x 10- 


4 


1.77 x 10 


-3 




6.0150 


3.03 x 10" 


5 




8.1464 x 10" 


4 


3.03 x 10 


-5 




6.0151 


2.18 x 1(T 


3 




8.1465 x 10- 


4 


2.77 x 10 


-3 



and P e < 1.10 x 10 5 . TableQ]shows how small variation of the sweep parameters alters Tr P. 
As with the other one-qubit gates, performance is most sensitive to variation in rj 4 , and the 



Table 4. Variation of Tr P for the NOT gate when the TRP sweep parameters are altered slightly 
from their best performance values. The columns to the left of center have ry 4 = 2.9277 x 10~ 4 
and those to the right have A = 7.3205. 



V4 


A 


TrP 


A 






TrP 




2.9277 x 10" 4 


7.3204 


1.12 x 10" b 


7.3205 


2.9276 x 10" 


-4 


1.23 x 10" 


3 




7.3205 


1.10 x 10~ 5 




2.9277 x 10" 


-4 


1.10 x 10- 


5 




7.3206 


1.22 x 10" 5 




2.9278 x 10" 


-4 


1.23 x 10" 


3 



sweep parameters must be controllable to 5 significant figures for the gate error probability 
P e to fall below the accuracy threshold P a ~ 10~ 4 (see Section |SJ). 

3.2 Modified Controlled- Phase Gate V CP 

We complete the TRP implementation of the universal set Q u by showing how the modified 
controlled-phase gate V CP can be produced. As shown in Appendix A, the two-qubit Hamil- 
tonian H 2 {t) used to implement V CP depends on two sets of parameters. The first set (A, jj 4 ) 
consists of the now familiar TRP sweep parameters, while the second set (^1,^2,^3,^4,04) 
consists of parameters such as coupling constants and frequency-related shifts and differences 
that appear in #2(7"). We stress that only the first set are directly related to the TRP sweeps. 
Unlike the situation encountered in Section [3.1i best gate performance does not require high 
precision control of the TRP sweep parameters A and rj A . Instead, the critical parameters for 
gate performance will turn out to be d\, di, and C4. As this is the first time we present results 
for V C p, we also include the unitary gate U a produced by our best-case choice of parameters. 

The target gate V CP has real and imaginary parts (see eq. (|3T)]0 



Re (Vo 



Im {V c 



(I 

1 



y 



\ 



-1 

1 / 



/ \ 





y ) 



(38) 
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The parameters 



A = 5.1 

7? 4 = 2.4 x 1CT 4 



d% = 11.702 
d 2 = -2.6 
d 3 = -0.41 
d 4 = 6.6650 



c 4 = 5.0003 



produce the gate U a with 



Re (Ua) 



Im(U a ) 



( 0.9998 
-0.0154 
0.0042 
\ -0.0026 

/ 0.0052 
-0.0109 
0.0030 
V -0.0018 



0.0155 
0.9997 
-0.0002 
-0.0021 

-0.0108 
0.0064 
0.0084 
0.0068 



0.0041 
-0.0003 
-0.9999 
-0.0037 

-0.0031 
-0.0084 
0.0060 
0.0079 



0.0028 \ 
0.0021 
-0.0038 
0.9999 / 

-0.0017 \ 

0.0068 
-0.0079 

0.0026 / 



(39) 



All two-qubit simulations used a dimensionless inversion time tq — 120.00. From U a and V CP 
it follows that TrP = 1.27 x 10~ 3 , the gate fidelity T Vcp = 0.9996 83, and P e < 1.27 x 10~ 3 . 
Table [5] shows how gate performance varies when either A or ?/ 4 is altered slightly. Notice 

Table 5. Variation of TrP for the modified controlled-phase gate Vcp when the TRP sweep 
parameters are altered slightly from their best performance values. The columns to the left (right) 
of center vary A (174), while holding all other parameters fixed. 



A 


TrP 




??4 




TrP 




5.0 


2.70 x 10- 


:s 


2.3 x 10- 


4 


1.46 x 10 


-3 


5.1 


1.27 x 10~ 


3 


2.4 x 10" 


-4 


1.27 x 10 


-3 


5.2 


2.10 x 10" 


3 


2.5 x 10" 


-4 


1.35 x 10 


-3 



that the TRP sweep parameters only need to be controlled to two significant figures. A 
similar situation occurs for the parameters di and cfa, though in the interests of brevity, 
we will not include Tables to show this. Table [6] shows how gate performance varies when 
either d\ or dn is varied slightly. We see that these parameters need to be controlled to five 

Table 6. Variation of TrP for the modified controlled-phase gate Vcp when the parameters d\ 
and d,4 are altered slightly from their best performance values. The columns to the left (right) of 
center vary d\ (di), while holding all other parameters fixed. 



di 


TrP 




G?4 


TrP 




11.699 


1.41 x 10" 


■> 


6.6647 


1.31 x 10- 


■1 


11.700 


7.63 x 10" 


3 


6.6648 


6.35 x 10- 


3 


11.701 


3.36 x 10- 


3 


6.6649 


2.40 x 10" 


3 


11.702 


1.27 x 10- 


3 


6.6650 


1.27 x 10" 


3 


11.703 


1.43 x 10" 


3 


6.6651 


2.97 x 10" 


3 


11.704 


3.79 x 10- 


3 


6.6652 


7.59 x 10- 


3 


11.705 


8.27 x 10" 


3 


6.6653 


1.50 x 10" 


2 



significant figures, although a small amount of uncertainty in the fifth significant figure will 
not drastically damage performance. Recall from Appendix A that d\ is the dimensionless 
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version of the difference in qubit Larmor frequencies, and is the dimcnsionless Ising coupling 
constant. Finally, Table [7] shows how Tr P varies when C4 changes slightly. We see that best 

Table 7. Variation of Tr P for the modified controlled-phase gate Vcp when the parameter C4 is 
varied slightly from its best performance value. The columns to the left (right) of center vary C4 
in the fifth (fourth) significant figure, while holding all other parameters fixed. 



c 4 


TrP 




C 4 


TrP 




5.0000 


1.98 x 10" 


-3 


4.999 


1.50 x 10" 


2 


5.0001 


1.55 x 10" 


-3 


5.000 


1.98 x 10" 


3 


5.0002 


1.36 x 10" 


-3 


5.001 


5.48 x 10" 


3 


5.0003 


1.27 x 10" 


-3 








5.0004 


1.38 x 10" 


-3 








5.0005 


1.65 x 10" 


-3 








5.0006 


2.11 x 10" 


-3 









performance is not seriously compromised if C4 has a small uncertainty in its fifth significant 
figure. Interestingly, if C4 can be controlled to four significant figures, we can come very close 
to best performance. As shown in Appendix A, C4 is a dimensionless degeneracy-breaking 
parameter. Although gate performance for V CP is slightly more robust than for the one-qubit 
gates of Section 13.11 we have not yet been able to find a combination of sweep parameters 
and two-qubit interaction that yields P e < 10~ 4 . See Section [5] for further discussion. 

4 Realizing TRP in Superconducting Qubit Systems 

Here we demonstrate how TRP sweeps can be applied to superconducting (SC) qubit systems 
[26l [27] . We first consider a SC charge qubit in Section [4.11 then go on to flux qubits in 
Section 14.21 Two realizations of a flux qubit are considered: first the rf-SQUID qubit in 
Section 14.2. 1[ then the persistent-current qubit in Section 14.2.21 For each of these qubit 
realizations, the demonstration proceeds in 3 steps: (i) the appropriate one-qubit Hamiltonian 
is introduced; (ii) the coefficients of a z and a x in this Hamiltonian are identified with at and 
b cos cf>t rp (t), respectively; and (iii) the rotating-wave approximation is invoked to obtain the 
one-qubit TRP Hamiltonian (see eqs. (p} and In each demonstration, the second step 

establishes the link between the theoretical parameters (a,b, (j} tr p(t)) and the experimental 
control fields acting on the SC qubit. 

4-1 Superconducting Charge Qubit 

Figure Q] shows a quantum circuit that can act as a charge qubit with adjustable Josephson 
coupling EH1 ED] - The circuit consists of a superconducting (SC) island that is connected 
to a SC electrode via a dc-SQUID that is threaded by an external magnetic flux & x , as well 
as to a gate voltage V g through a capacitor C g . The dc-SQUID is designed to have low self- 
inductance and is made up of two identical Josephson junctions (JJ), each with Josephson 
coupling energy Ej and capacitance Cj. The SC energy gap A is assumed to be the largest 
energy scale for the circuit dynamics so that at sufficiently low temperature only Cooper pairs 
can tunnel through the JJs. Denoting the number operator for Cooper pairs on the island by 
n and the phase of the SC order parameter by 8, the Hamiltonian for the SC island is 

H = AE c {n -n g f - Ej{<5> x )cos6. (40) 
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E°j,Cj 






SC Island 




r 1/9 



Fig. 1. Quantum circuit for a charge qubit with tunable Josephson coupling. The superconducting 
(SC) island is connected to a SC electrode through a dc-SQUID that is threaded by an external 
magnetic flux <& x . Each cross X in the dc-SQUID represents a Josephson junction. The island is 
also coupled to a gate voltage V g through a capacitor C g . 



Here: (i) n g — C g V g /2e; (ii) E c = 0.5e 2 /(C g + 2Cj) is the charging energy; and (hi) the flux- 
controlled Josephson coupling energy is Ej($ x ) = 2Ej cos (7r$ x /$o)> where $0 = h/2e is the 
flux quantum. The number operator h and the phase 9 are canonically conjugate variables 
subject to the uncertainty relation AnAd > 1 [31 . In the charging limit [E c Ej) H is 
dominated by the charging term in eq. (|40[) . In this case the eigenstates \n) of n form a 
convenient basis with which to represent H. Because of Cooper pair tunneling, the Josephson 
coupling term in H will have matrix elements connecting the states \n) <-> \n + 1). It follows 
from these remarks that H can be written as 

H = J2{^E c (n~n g ) 2 \n)(n\~^^l[\n)(n + l\ + \n + l){n\}\. (41) 

Note that in the absence of the Josephson coupling term, and for n g = 1/2, the states |0) and 
|1) have degenerate energies: Eq = E\ = E c . In fact, a level-crossing occurs at this value of n g . 
The presence of the Josephson coupling in H causes an avoided crossing to occur at n g = 1/2, 
opening an energy gap Ai? i = Ej($ x ) <C E c . Near the degeneracy point (n g — 1/2), 
these two states have the smallest energy eigenvalues. Thus at low temperature, and for 
low frequency gate voltage V g (t) and external magnetic flux <fr x (t), the circuit dynamics near 
the degeneracy point is restricted to the subspace spanned by the states |0) and |1). Under 
these conditions, the only terms in H that are dynamically relevant are those that act on 
this subspace. Truncating H so that only those terms are kept gives the SC charge qubit 
Hamiltonian 

H oq = -la-B, (42) 

where 

B z = 4E c {l-2n g ) 
B x = Ej($ x ) 

By = 0. (43) 
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This completes the first step of the demonstration. 

To apply a TRP sweep to a charge qubit we must require that 



B z 

~2 = 

-y = b cos (ptrp, (44) 

where (ptrp = (2/n)Bt n is the TRP twist profile. This is step 2 in the demonstration. With 
these assignments, eq. (|42j) becomes 

H cq = —ata z — b cos 4> trp a x . (45) 

In the rotating wave approximation (step 3) this becomes 

H cq =-ata z -b cos cj) trp a x ~b sin cj) trp a y . (46) 

Comparing this with eqs. (p} and (|2|), we see that eq. (|46l) is the Hamiltonian for a qubit 
interacting with a TRP sweep. Plugging eqs. ([4*3"j) into eqs. (j4~4"|) , and recalling that Ej(& x ) = 
2E ( jcos(7r^ x /^o) gives 

= (^j (f>tr P (t) (48) 
£° = b. (49) 

Eqs. (|47|) and ([48]) specify the time dependence that V g (t) and "^(i) must have, respectively, 
for a TRP sweep to be applied to a SC charge qubit, and eq. (|49|) links the sweep parameter b 
to the Josephson coupling energy Ej. The TRP sweep parameters a, B, and To can then be 
found from eq. (f32|) for given values of A, ?y 4 , and tq. Varying the gate voltage V g (t) according 
to eq. (|47|) produces the inversion of the TRP control field F(£) in eq. ((2]), while varying 
the flux <& x (t) through the dc-SQUID according to eq. (j48|) . together with eq. (j49|) and the 
rotating-wave approximation, produces its twisting in the x-y plane. 

4-2 Superconducting Flux Qubit 

We examine two flux qubit proposals: (i) the rf-SQUID qubit (Section 14.2. and (ii) the 
persistent-current qubit (Section 14.2.21) . 

4.2.1 rf-SQUID Qubit 

In an rf-SQUID a single Josephson junction (JJ) interrupts a superconducting (SC) loop 
that is threaded by an external magnetic flux Q x . Because the loop has non- vanishing self- 
inductance L, a secondary flux $ s = Li s is present whenever a supercurrent i s circulates 
around the loop. The total flux $ = & x + $ s determines: (i) the phase difference ip across 
the JJ via ip — 27r$/<I > o (mod 27r); and (ii) the supercurrent i s via i s = i c sin(27r<i>/(f>o)- Here 
i c is the critical current for the junction, and <i>o = h/2e is the flux quantum. The potential 
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energy for the rf-SQUID is the sum of the Josephson coupling energy and the magnetic energy 
associated with the secondary flux <E> S = $ — <fr x : 



here Ej is the coupling energy. £/($) has a number of important properties. First, U(§>) 
forms a double-well potential for <&-values near &o/2 when: (i) Q x <&o/2; and (ii) the self- 
inductance L is large enough to cause (3l = E® / / An 2 L) > 1. The minima occur at 4>± = 
$o/2±(5$±, and the supercurrents corresponding to these minima, i± = =Fi c sin(27rJ$±/$o)i 
have opposite circulations. When & x = $0/2, the double-well potential is symmetric; the 
central barrier has a maximum at $ = $o/2, and the barrier height at maximum is Ej. 
When $> x = $o/2 + 8& x , the double- well potential is asymmetric, and for small S$ x , the 
barrier height at maximum remains of order Ej. 

The quantum degree of freedom for an rf-SQUID is the total flux or equivalently, the 
supercurrent i s . In the absence of tunneling through the central barrier, the groundstate is 
doubly degenerate and the energy eigenstates \L) and \R) are localized, respectively, about 
the left and right minima of the double-well potential. Tunneling splits the degeneracy and 
the new energy eigenstates \E a ) and \E a ) are, respectively, symmetric and antisymetric linear 
combinations of \L) and \R). At sufficiently low temperature the rf-SQUID is limited to the 
subspace of states spanned by \E S ) and \E a ). Bocko et al. [32] were the first to propose using 
a low-temperature rf-SQUID as a qubit. The computational basis states (CBS) |0) and |1) 
are identified with the states \R) and \L), respectively, and are defined to be eigenstates of the 
operator a z : cr z \i) — (— with i = 0, 1. For $> x w $o/2, the asymmetry of the double-well 
potential introduces a bias e in the minima of the double-well: 



In the absence of tunneling, the energy E (Ex)) of the state \R) = |0) (\L) = |1)) is well 
approximated by the sum of [/($+) (£/($_)) and the groundstate energy of a harmonic 
oscillator with frequency determined by the double- well curvature at $ = ($ = $-). For 
small asymmetry, the energy difference E\ — Eq for the CBS is then e. To account for this, 
the flux qubit Hamiltonian Hf\ includes the term — (e/2) <7 Z , where the zero of energy has 
been set at {Eq + E\)/2. Note that this term in Hfi can be varied by varying the flux & x 
which alters the asymmetry of the double-well, and thus the bias e. Tunneling between wells 
causes a bit-flip operation |0) <-> |1) to be applied to the qubit state. If [fit, — 1) <C 1, eq. ([ST]) 
indicates that |e| <§C Ej, and so the energies =F e /2 of the two lowest energy eigenstates are 
well below the barrier maximum. WKB theory can thus be used to determine the amplitude 
ti r /2 for the qubit to tunnel from the left to the right well (the factor of 1/2 is introduced for 
convenience). On the basis of the above remarks, we are led to the flux qubit Hamiltonian 




(50) 




(51) 




(52) 
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where 



B z = e 

B x = tl r 

B y = 0. (53) 

The tunneling amplitude U r /2 is highly sensitive to the barrier height which, as noted earlier, 
is of order Ej. If the JJ in the rf-SQUID is replaced by a dc-SQUID threaded by an exter- 
nal flux $> x (as in Section l4~Tj) . the coupling energy becomes adjustable, Ej — > Ej(Q x ) = 
2E°j cos(ir<b x /Q ), and $ x can be used to vary B x — ti r . With this modification, both terms 
in Hfi are experimentally controllable. This leads us to the quantum circuit in Figure [2] for 
an rf-SQUID flux qubit. The JJs inthe dc-SQUID are identical and have coupling energy Ej 



E°j,C 




E°j,C 



Fig. 2. Quantum circuit for an rf-SQUID flux qubit with tunable Josephson coupling. The larger 
superconducting loop is threaded by an external magnetic flux <fr x , while the smaller loop in the dc- 
SQUID is threaded by the external flux & x . Each cross X in the dc-SQUID represents a Josephson 
junction with coupling energy Ej and capacitance C. 



and capacitance C. The WKB expression for ti r /2 is [33] 



tlr TllQ* 

~2 ~ ~2^ 



exp 



df y/2M nn \E-U\ 



(54) 



Here is the attempt frequency; / = 2tt [($/$q) — (1/2)]; M nn = 2C; and E is the energy. 
Following Ref. [33] we find that 



(h 



1 



LC ' 

and denoting the argument of the exponential in eq. (|54j) by /, we find that 



I = 



ILC 



((3 L -lf 2 Ej[<S> x 



Eqs. ([S"2"]) - (p)()|) complete the first step of our demonstration. 

To produce a TRP sweep the time-dependence of <fr x and <t x must be such that 



b cost 



at = 



trp 



B_z_ 
2 

B x 
~2 



(55) 



(56) 



(57) 
(58) 
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In principle, this completes step 2 of the demonstration. In actuality, a bit more work is needed 
to fully establish the link between the theory parameters and the control fluxes appearing in 
Figure [21 To that end, substituting for B z using eqs. (|5"Tj) and ([53"]) in eq. (|5T|) gives 



at = 2tt 







I $0 


4} 



Ej($ x )y/6(p L -l) 



(59) 



where E° -> £j($ x ) due to the dc-SQUID in Figured and (3 L = Ej($ x )/ (<Z> 2 /4ir 2 L) . It 
will be seen below that eq. (|58p requires <& x to be time-dependent, and yet we would like the 
term in the square bracket in eq. (|59[) to be constant. To achieve this we write 



(60) 



where $° is time-independent and will be specified below, and we require that n5$ x /$o <C 1 
Then we can write 

Ej(* s ) = Ej(*°) + SEj 



and 



where 



Ej(¥ x )=2E jC o S 



$0 



SEj = -£j($°)tan 



$ 



$0 



and 



01 = 



\4tt 2 L J 



S(3 L 



SEj 



\4:Tr 2 L J 



(61) 
(62) 

(63) 
(64) 



We now require that \SEj\ -C Ej($ x ) and 5(3l <C /3° — 1 so that eq. ([59")) can be written as 

(65) 







I *0 


4} 



Solving for <f> x (i) gives 



*.(«) = £ 



at 



(66) 



By requiring that the flux ^(i) satisfy eq. (|66|) . we insure that the a z term in -ffyi has the 
appropriate at coefficient needed for a TRP sweep. Thus appropriate variation of the flux 
$ x (t) through the primary loop in Figure [2] produces the inversion of the z-component of the 
control field F(t) in eq. @. Next we use eq. I|58p to determine the flux $ x (t) through the 
dc-SQUID in Figure [2] Using eqs. (|53l ) -([56 |) in eq. (58|) gives 



b cos< 



'trp 



2tt 



■ w» e 



(67) 
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Writing I = I + 61, where 61 -C I and 



p(/3°-l) 3/2 ^($0) 



SI = In 



3 6j3 L 



5Ej 



2(/3°-l) Ej^l) 



allows eq. (j6"7| to be written as 

h 



b cos ( 



Ztt 



l-Io 



3 6f3 L 



SE, 



2 (Pi 1) Ej{¥ x ) 
Note that cj» = uj® + 5uj* has been written as w° in eq. (l70|) , where 



(68) 
(69) 

(70) 



(71) 



The reason for this is that the 61 contribution to eq. (|70[) dominates the 5uj* contribution, 
and so the latter contribution can be safely discarded. Using eqs. (J63J) and (|64|) in eq. ([70|) . 
along with some algebra, gives 



( b cos 4> trp ) 



2tt 



Jo (5/3° -2) SEj 
2 (/?£-!) 



(72) 



Finally, using eq. ([63]) again, and solving for n6<P x /§ gives 

ir6$ x 



= C cos < 



(73) 



where 



C 



D 




CSC 



$0 



(5/30-2)7^1 



exp [ 7 ] 



LCj (5(31-2)^1^1 \E°j 



(74) 



Typical values for the parameters in eqs. (|74[) are: — 1 = 0.1; \JLC = 1 ns; and Ej/h = 
100 GHz. With these values, and choosing b/fi = 400Hz, and 

$°/$o = 1/2 - c/tt, (75) 

with e = 0.25 gives 



0.32 GHz 
12.7. 
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From these values we find that 



C = 
D = 



2.9 x 1(T 3 
1.1 x 1(T 3 , 



and from eq. l[73|) . it follows that ir5$ x /<S>o w 10 -3 « 1 as required. The TRP sweep 
parameters a, B, and T are then determined from eq. (|32p for given values of A, ry 4 , and r . 
Thus, with <I> X given by eq. ([66]) and <& x given by eqs. ([60]) . (|73|) . and ([75]) . we arrive at the 
flux qubit Hamiltonian is 

Hfi = —ata z — b cos 4>trp<?x- 

We see that <fr x (t) produces the inversion in the TRP control field F(i) in eq. J5J), while $ x (*) 
will be seen momentarily to give rise to its twisting in the x-y plane. This finally completes 
step 2. In the rotating wave approximation (step 3), Hfi becomes 



Hfi = —at <t z — b cos (f>trp &x — b sin <j) trp a y , 
which is the Hamiltonian for a qubit acted on by a TRP sweep (eqs. ((I) and (J2|)). 



(76) 



4-2.2 Persistent- Current Qubit 

The focus here is the 4-junction persistent-current qubit introduced in Refs. [33l [34l [35] . The 
quantum circuit for this qubit is shown in Figure [3l It consists of a loop with negligible 



Va 



V b 



CgA, VgA 



SCI 



Eji,Ci :c 




Ej3 , C 3 



CgB,VgB 



SC2 



X E J2 ,C 2 



Fig. 3. Quantum circuit for a 4-Josephson junction persistent-current qubit. Junctions 1 and 2 are 
identical and have coupling energy Eji = Ej2 = E°j and capacitance C\ = Ci = C. Junctions 3 
and 4 are also identical, with coupling energy Ej% = Ej^ = /3-Ej and capacitance C3 = C4 = j3C. 
Superconducting islands SCI and SC2 are, respectively, connected to gate voltages Va and Vb 
through gate capacitors C g A and C g B, where C g \ = C g B = jC. The potential difference across 
each gate capacitor is V g i, where I = A, B. Finally, the dc-SQUID loop is threaded by a flux $2, 
while the lower loop is threaded by a flux <J>i. 



self- inductance that is threaded by an external magnetic flux $1. Because of the small self- 
inductance, a supercurrent i s circulating around the loop produces negligible secondary flux 
and so the total flux through the loop <& is given by the external flux $1. The loop is 
interrupted by two identical Josephson junctions 1 and 2; a dc-SQUID containing two identical 
junctions 3 and 4 whose loop is threaded by a flux $2; and two superconducting islands 
SCI and SC2. Junctions 1 and 2 both have coupling energy Ej and capacitance C, while 
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junctions 3 and 4 have coupling energy (3Ej and capacitance (3C. Superconducting island SCI 
(SC2) is connected to a gate voltage Va (Vb) through a gate capacitor C 9 a {C 9 b)- The gate 
capacitors are identical with C g A = C g s — lC 1 and capacitor C g j has a potential difference 
V g i across it, where I = A,B. The circuit is assumed to be at sufficiently low temperature 
that only supercurrents flow in it. 

It is conventional to introduce the magnetic frustration /j = <f>i/$0: where $o = h/2e is 
the flux quantum and i = 1, 2. Denoting the Josephson phase difference across the j th junction 
by ifj, fluxoid quantization around the dc-SQUID loop in Figure [3] requires (p^ — ips = — 2nf2, 
and around the lower loop requires (fx — <p% +if3 — — 2itf\. Adding up the Josephson coupling 
energy for each junction and using the fluxoid quantization relations gives the total Josephson 
energy U : 

-^o = 2 + 2/3 - 2 cos (p p cos tp m - {2(3 cos (7r/ )} cos (2tt/ 6 + 2tp m ) , (77) 

where (p p = (cpi + <p 2 ) /2; ip m = (ipi - ip 2 ) /2; f a = / 2 ; and f b = fa + (/ 2 /2). The Josephson 
energy £/ acts as the potential energy for the persistent-current qubit. It is periodic in /& 
with period 1, and is symmetric about fb = 1/2. At the degeneracy point /& = 1/2, a = 
2/3cos(7r/ a ) = 1/2, U has degenerate minima at (<p p ,(p m ) = (0,±(/?^), where cosip* n = l/(2a). 
Because C/ is periodic in ip p and ip m , the degenerate minima trace out a lattice whose primitive 
unit cell contains only one pair of degenerate minima. When fb is near 1/2, the pair of minima 
inside a unit cell are no longer degenerate. Within a unit cell, the potential U takes the form 
of a double well potential which is symmetric when /j, = 1/2, and is asymmetric when fb is 
near 1/2. Variation of a allows the potential energy landscape to be further modified. In 
particular, a can be used to suppress tunneling between minima located in different unit cells. 

In the quantum limit: (i) the fluxes tp p and (p m become quantum degrees of freedom with 
conjugate momenta P p — —ifid / dip p and P m = —i%d/d<p m \ and (ii) the Hamiltonian driving 
the quantum dynamics is 

p2 p2 

2M P 2M m J 

Here M p = ($ /2tt) 2 (1 + 7) (2C) and M m = ($o/2vr) 2 (1 + 4/3 + 7) (2(7). The energy levels 
form bands due to the periodicity of U, and the bands are symmetric about fb = 1/2. At 
sufficiently low temperature, and for fb near 1/2, the circuit is effectively restricted to the 
subspace spanned by the two lowest energy eigenstates \Eq) and \Ek). Persistent supercurrents 
i and i\ flow in \E ) and \Ei), respectively, and i = — ij. These two eigenstates constitute 
the CBS |0) and |1), and the bit values are encoded into the direction of circulation of the 
persistent-currents. In the tight-binding approximation the states \Eq) and \E±) are localized 
near the minima of U, and (as noted earlier) within a unit cell, U has the form of a double 
well potential. Focusing on the unit cell containing the minima (0, ±</3^J, and denoting 
the state localized about the minimum (0,+ip* n ) ((Q, — by |+) (|— )), inside this unit 
cell we have \E ) = [|+> + |-)]/v^ and = [|+) - |-)]/\/2. When a = 2f3cos{nf a ) is 
chosen appropriately, tunneling can be restricted to the two minima within a unit cell. In the 
absence of tunneling, |+) (|— )) is the groundstate of a qubit localized in the well centered 
about (0, +<Pm) ((0, — Altering fb alters the energies E± of the states |±), and so 

alters the bias F = (E + — E_)J2, Altering f a alters the barrier height, and so alters the 
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tunneling amplitude t. Thus, in the representation spanned by |±), the persistent-current 
qubit Hamiltonian H$i is 

' —F -t 
-t F 



H f2 = 



Since the CBS are the eigenstates of -ff/2, we can also write 



H 



J-2 



■VF 2 + t 2 a z , 



where a z — |0)(0| — Following Ref. [33], the quantum circuit is to operate near the 

degeneracy point at which /* = /| = 1/3. The operating point is chosen to have f° = f* + e\ 
and /J 1 = + £2- It can be shown that, at the operating point, the bias is F a = r±ei + r^ti 
and the tunneling amplitude is t a = t* + S2£2, where is the tunneling amplitude at the 
degeneracy point. Explicit formulas for n, T2, and s 2 appear in Ref. [53], though they will 
not be needed here. By varying /1 = /° + <5i and / 2 = /I + £2 about the operating point with 
time-dependent £1 and (5 2 , a unitary transformation can be applied to the qubit. Defining 
tan# = —t /F , the Hamiltonian Hf2 becomes 



H f2 = [-Z + Zi <5i + Z 2 6 2 } <J Z - [XkJi + X 2 <J 2 ] a,, 



where 



^0 = ^+^! 

Zi = n cos # 

Z 2 = ^ cos 9 — s 2 sin 9 Q 
Inserting typical parameter values, Hf2 becomes 



X\ = n sin 9 
X 2 = ^ sin o + s 2 cos 9 a 



= [ -0.025 + 4.0 5i + 2.1 62] o z 



\ 0.46 <5i + 0.41 <5 2 ] ct s 



(78) 



(79) 



(80) 



This completes step 1 of the demonstration. 

To produce a TRP sweep, we must require (step 2) that 

Zq — Z\8\ — Z 2 S 2 = at 

Xi Si + X 2 S 2 = b cos 4> trp - 



(81) 



Recalling that t — (b/a)r, and defining Z% — Ej Zi and Xi — EjXi, eq. (|5Tj) can be rewritten 



as 



z\ Si + z 2 S 2 

Xi 81 + X2 82 



b 



Solving for 81 and 82 gives 

81 = 

82 = 



X2 






G 




)- 


Xi 


zo - | 


' b ^ 


G 





-'0 



Z2 (b_ 

G U.°/ 

Zl (b_ 

G U.°/ 



^trp 



COS 4>trp, 



(82) 



22 



where G = X\Zi — xiz\. Typical values for Xi and z$ for i = 1, 2 can be read off from eq. (|80[) . 
Inserting these values into eq. (|82|) gives 



5 x ( r ) = 0.61 
<5 2 (t) = 0.68 



_6_ 



T - Z 



>!J ( jo ) ( ' <IS ""7' 



5.9 ( ) cos (j) trp - (83) 



Recall that a TRP sweep runs over times — to/2 < t < to/2, and that our simulations used 
t /2 = 40,60. In Ref. [33], Si ~ 10" 3 and 5 2 ~ 10~ 4 . To produce Si and <5 2 of this size 
using eq. (Ml, we require that b/E 7 ~ 10~ 4 and z ~ 10" 4 . With E Q j/h = 100 GHz, this 
requires b/h, Zo/fi ~ 10 MHz. The TRP sweep parameters a, B, and To are then found 
from eq. (|3"2"|) for given values of A, 774, and tq. Thus, with Si and S2 given by eqs. (|53"|) . the 
persistent-current qubit Hamiltonian becomes 

if/2 = — 0.t <J Z — b COS (f>trp Px- 

This completes step 2. In the rotating wave approximation (step 3), Hf2 becomes 

Hf 2 = —at a z — b cos 4> tr p &x — b sin 4> trp a y (84) 

which is the Hamiltonian for a qubit subjected to a TRP sweep (eqs. (HJ and @). We see 
that requiring the control fluxes Si(t) and S%{t) to satisfy eqs. ([55)1 causes a TRP sweep to be 
applied to the persistent-current qubit. Unlike with the two previous SC qubit realizations, we 
see from eq. (|5T1) that the inversion and twisting of the TRP control field F(t) is determined 
by linear combinations of the control fluxes Si and 82 ■ In the case of the persistent-current 
qubit, we cannot (in general) attribute inversion to one control field and twisting to the other. 

5 Discussion 

In this paper we have presented simulation results which suggest that TRP sweeps should be 
capable of implementing a universal set of quantum gates Q u that operate non-adiabatically 
and with high-fidelity. The one-qubit gates in Q u were seen to operate with gate error 
probabilities satisfying P e < 10~ 4 , and the two-qubit modified controllcd-phase gate with 
P e < 1.27 x 10~ 3 . Using the rough-and-ready estimate for the accuracy threshold, P a ~ 10 -4 , 
we see that: (i) the TRP-generated one-qubit gates have error probabilities that fall below 
this threshold value; and (ii) the two-qubit gate comes within an order-of-magnitude of it. 
As was noted in Section IU this high level of gate performance requires that the TRP sweep 
parameters be controllable to high precision. Finding a way to make TRP gates more robust 
is a major challenge that must be overcome if these sweeps are to become a viable means 
for universal control of a quantum computer. To that end, we are currently exploring the 
consequences of interlacing TRP sweeps with dynamical decoupling pulses [351 133 13H] to pro- 
duce an effective dynamics that preserves a group of symmetries G of the target gate U t . By 
removing the part of the TRP dynamics that does not commute with G, it is hoped that 
the resulting effective dynamics will yield gate performance that varies more slowly with the 
sweep parameters, and makes it easier for the minimization algorithms to carry out their 
work. 
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TRP has been experimentally realized in NMR systems [TSJ, [T7] , and these systems have 
allowed the controllable quantum interference effects that were predicted to arise from multiple 
passes through resonance per TRP sweep [TB] to be observed [17] , A demonstration of how 
TRP sweeps can be applied to atomic systems using electric fields has also been given [18] . In 
this paper we have shown how TRP can be applied to both superconducting charge and flux 
qubits. In such superconducting qubit systems, a qubit is realized by a mesoscopic circuit 
and the control fields are gate voltages and highly localized magnetic fields that are applied 
directly to the circuit. These systems thus allow a TRP sweep to be applied to an individual 
qubit, unlike in NMR systems where the sweeps are applied to an ensemble of qubits. In 
the NMR setting, variation of the rf-magnetic field amplitude over the ensemble causes the 
largest problem for high precision control of the TRP sweep parameters [TB]. Note that this 
difficulty does not arise with superconducting qubits. It would be interesting to work out an 
adaptation of the NMR experiment that observed the TRP quantum interference effects [T7] 
to a superconducting qubit system. Since these interference effects arc a direct consequence 
of the temporal phase coherence of the qubit wave function, the adapted experiment would 
provide a new, independent demonstration of quantum coherence in superconducting qubit 
systems. Such a translation of the NMR experiment is currently underway. Ref. jT8] described 
a state tomography experiment that would test the TRP gate simulation results by measuring 
the output density matrix p exp — U a \tpo) (tpo\Ul resulting from an initial state |^>o)j f° r each 
of the TRP generated gates U a presented in Section [3] Associated with each sweep is a target 
gate Ut and a corresponding target density matrix p t = Ut\ipo) {ipo\Ul . Having measured p exp , 
the gate fidelity could be calculated and compared with the fidelity obtained from the TRP 
gate simulations. There are now three possible physical systems where such an experiment 
could be carried out. For alternative applications of temporal phase coherence and rapid 
passage to quantum computing, see Refs. [39] and [lQl|4T], respectively. 
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Appendix A 

Here we derive the dimensionless Hamiltonian H2(t) introduced in Section 12.21 It proves 
convenient to adopt the language of NMR, although a more general discussion is possible. 

Consider a two-qubit system in which each qubit is Zeeman-coupled to an external mag- 
netic field B(i), and the two qubits interact through the Ising interaction. As noted in Sec- 
tion [^21 the Ising interaction was chosen because of its simplicity, and because it is present 
in many physical systems. It is a simple matter to alter the following arguments to include a 
different two-qubit interaction. Our starting point is thus the Hamiltonian 



where ji is the gyromagnetic ratio for qubit i, and J is the Ising interaction coupling con- 
stant. In the lab frame, B(t) has a static component BqZ and a time-varying component 
2B r f cos 4> r f(t) x. In the rotating wave approximation B(i) reduces to 



Introducing = ^Bq and w\ = jiB r f (i = 1, 2), and inserting eq. (|A.2|) into eq. (jA.lj) gives 




(A.l) 



B (f ) = BqZ + B r f cos (f> r f (i) x — B r f sin cf> r f (t)y. 



(A.2) 
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h 



E 




{cos 4> r fO-ix - sin <j>rfO- iy } ~-Jo- lz a 2z - (A.3) 



i=l 



Transformation to the detector frame is done via the unitary operator 



U{t) = exp [ (# det (t)/2) (a lz + a 2z ) ) . 



The Hamiltonian in the detector frame is then [19] 




^ (- y + <t>detj o- lz )— {cos (<f) det - cf> r f) a ix + sin (<j) det - <t> rf ) a iy } 

i=i L 

-7T J (J lz< J 2z- 



(A-4) 
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As explained in Section 12.11 to produce a TRP sweep in the detector frame it is necessary to 
sweep <f>det and r / through a Larmor resonance frequency. We choose (somewhat arbitrarily) 
to sweep through the Larmor frequency oj 2 : 

■ 2at . 

<Pdet = uj 2 + — + A 

h 

(j) r f = <fi det - 4 . (A. 5) 



Here 04 (i) = (l/2)Bt 4 is the twist profile for quartic TRP, and we have introduced a frequency 
shift parameter A whose value will be determined by the sweep parameter optimization 
procedure of Section [231 Inserting eqs. (IA.5|) into eq. (IA.4I) . and introducing 5uj = w\ — ui 2 
and bi = huj r J /2 (i — 1, 2), we find 
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(A.6) 



We see that both qubits are acted on by a quartic TRP sweep in the detector frame. In keeping 
with our earlier choice of sweeping through the Larmor resonance of the second qubit, we use 
b 2 in the definitions of the dimensionless time r, inversion rate A, and twist strength 7/ 4 : 



fiB 
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(A.7) 
(A.8) 
(A.9) 



Since Yi 2 {t)/h has units of inverse-time, and b 2 /a has units of time (eq. (|A.7jl ). multiplying 
eq. (|A.6|) by b 2 /a and using eqs. (|A.7|) — (|A.9[) gives the dimensionless two-qubit Hamiltonian 
H 2 (r): 



H 2 (r) 
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where 




(A.ll) 



As noted in Section 12.21 H2 (t) has a degeneracy in the resonance frequency of the energy 
level pairs (E± «-> E%) and (£3 <-> E4). To break this degeneracy we add the term 

AH = c i \E i (T))(E i ( T )\ (A.12) 

to 62(7"), where \E±(t)) is the instantaneous energy eigenstate of 62(1") with eigenvalue £4(7"). 
Our final Hamiltonian is then 

fla(r) = H 2 (t) + AH (A. 13) 

which is the Hamiltonian given in eq. (|19[) . We see that i^C?") depends on the TRP sweep 
parameters (A, 7? 4 ), as well as on the parameters (di, . . . , (^4) and C4. From eq. (|A.11[) we see 
that di, d 2 , ds, and d4 are the dimensionless versions of, respectively, the Larmor frequency 
difference Sto = lu± — o; 2 , the frequency shift parameter A, the ratio 61/62 = 71/72, and the 
Ising coupling constant J. 
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